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We present the numerical results for the dynamics of a single chain in steady shear flow. The 
chain is represented by a bead-spring model, and the smoothed profile method is used to accurately 
account for the effects of thermal fluctuations and hydrodynamic interactions acting on beads due to 
host fluids. It was observed that the chain undergoes tumbling motions and that its dimensionless 
frequency F = Gnr^a'^ u / ksT depends only on the Peclet number Pe with a power law F oc Pe", 
where fee is the Boltzmann constant, T is the temperature, and a is the diameter of the beads. The 
exponent a clearly changes from 2/3 to 1 around the critical Peclet number. Pec, indicating that the 
crossover reflects the competition of thermal fluctuation and shear flow. The presented numerical 
results agree well with our theoretical analysis based on Jeffrey's work. 

PACS numbers: 83.80.Rs, 47.57.Ng, 82.20.Wt 



I. INTRODUCTION 

The dynamics of solid particles dispersed in host fluids 
is an important problem in many different fields of science 
and engineering. The macroscopic properties of such dis- 
persions (the elastic modulus, viscosity, and thermal and 
electric conductivities) greatly depend on the dynamics 
of the particles in the host fluids. In equilibrium states, 
the dynamics of small dispersed particles are strongly af- 
fected by the thermal fluctuations of their host fluids. 
When flow is imposed, the dynamics are also affected 
by the How of the host fluids. Because it is difficult to 
experimentally analyze these complex particle dynamics, 
which are coupled both to thermal fluctuations and to 
fluid flow, numerical simulations are particularly impor- 
tant for understanding the properties of particle disper- 
sions in detail 

For a single Brownian chain fluctuating in the shear 
flow of a Newtonian fluid, it has been suggested that the 
tumbling frequency v is proportional to the shear rate 
7^/^ [2nyl- This has been experimentally confirmed by 
Schroeder, Teixeira, Shaqfeh, and Chu for the dy- 
namics of individual DNA molecules in a linear shear 
flow. Primarily, Smith, Babcock, and Chu [loj mea- 
sured the power spectral density (PSD) and the probabil- 
ity distribution function (PDF) of the extension length 
of each DNA molecule for various Weissenberg numbers 
Wi = 7r, where 7 is the shear rate and r is the relax- 
ation time of the chain orientation. The PSD of polymer 
extension exhibits no peaks. In subsequent experiments 
0, Q, however, the focus has been put on the PSD of the 
orientation angle ip, where ip = when the DNA molecule 
lies perfectly in the flow direction. These experimental 
results support a simple power law, i/r oc Wi^^^ , where v 
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is the peak frequency of the PSD. The relaxation time t 
is considered to be a constant if the temperature is con- 
stant. This leads to i/ oc 7^/'^. Similar results have also 
been obtained in other experiments [1, Q jiiumerical sim- 
ulations [H, Q , and theoretical analyses ■ Although 
the DNA molecules mentioned above can be considered 
flexible chains, a quite similar power law was obtained 
using numerical simulations for a Brownian linear rigid 
rod as well 

For a single non-Brownian (thermally non-fluctuating) 
flexible chain in shear flow, the tumbling frequency v is 
expected to be proportional to the shear rate 7 A 
single non-Brownian rigid rod with a flnite aspect ratio is 
known to exhibit a cyclic tumbling motion in shear flow, 
as described by Jeffrey's equation, v oc j 

From the above experimental findings, one would ex- 
pect to observe a crossover from Brownian (j/ oc 7^^'^) 
to non-Brownian (z/ oc 7) behavior with increasing shear 
rate; however, such a clear crossover has not yet been re- 
ported. This crossover has not yet been successfully pre- 
dicted by previous numerical [sl, Q and theoretical [3-@ 
studies where a dispersed chain is treated as an end-to- 
end vector, namely, as an infinitely thin line. Therefore, 
the rotational motion of the chain cannot be sustained 
across cp = without thermal fluctuations. When a ther- 
mal fluctuation exists, the orientation of the thin line can 
fluctuate around ip = 0. This leads to tumbling motions 
even for a thin line; however, the frequency of the tum- 
bling motion v is always proportional to 7^/"^ , regardless 
of the shear rate 0. 

Wc thus aim to analyze this crossover by use of a direct 
numerical simulation (DNS) approach. In the present 
study, wc simulated the tumbling motion of a chain using 
a smoothed profile method (SPM) that accurately takes 
into account thermal fluctuations and hydrodynamic in- 
teractions The chain is represented by a bead- 
spring model, wherein each bead is modeled as a spheri- 
cal object with a finite radius a and undergoes free rota- 
tion. Rigid rods or fiexiblc chains are represented with or 
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without a constraint force on bond bending. Apart from 
the previous numerical models, the presented bead-spring 
model naturally takes into account the finite thickness 
of the experimentally used chains or rods. A theoreti- 
cal analysis has also been developed to understand the 
mechanisms underlying the crossover. 



II. METHODS 

A. MODEL 

We solve the dynamics of a single chain in a Newtonian 
solvent using SPM [isl - fTsj . In this method, boundaries 
between solid particles and solvents are replaced with 
a continuous interface by assuming a smoothed profile. 
This enables us to calculate hydrodynamic interactions 
both efficiently and accurately, without neglecting many- 
body interactions. The equation governing a solvent with 
a density pf and a shear viscosity 77 is a modified Navier- 
Stokes equation. 



{ dv , ^ - , ^ 



-VP + VV V + Pi4>fp + /shear (1) 



with the incompressible condition \J ' — 0? where 
^(r, t) and p(r, t) are the velocity and pressure fields of 
the solvent, respectively. A smoothed profile function 
< distinguishes between the fluid and par- 

ticle domains, yielding (/) = 1 in the particle domain and 
= in the fluid domain. These domains are separated 
by thin interstitial regions, the thicknesses of which are 
characterized by The body force 4>fv introduced to 
ensure the rigidity of the particles and the appropriate 
non-slip boundary condition at the fluid/particle inter- 
face. The mathematical expressions for (j) and (pfp are 
detailed in previous papers [l^ [l^ . The external force 
.^hear IS introduccd to maintain a linear shear [T9| . ex- 
pressed by 



-iV (0 < 2/ < Ly/A), 

-7(-y+^) {Ly/A<y<iLy/A), (2) 
-i{y-Ly) i3Ly/4<y<Ly), 

where 7 is the shear rate and Ly is the system size in the 
y-direction. 

In the present study, the chain is represented as either 
a rigid rod or a flexible chain. We use a bead-spring 
model consisting of N beads in a single chain. The bead 
size is sufficient to fit several mesh units. Therefore, it 
is necessary to consider the torque exerted on the bead. 
The motion of the ith bead is governed by the following 
Newton's and Euler's equations of motion with stochastic 
forces: 



d - 
at 



dt 



dt 



(3) 



(4) 



where Ri, Vi and Qi are the position, translational ve- 
locity, and rotational velocity of the beads, respectively. 
Mi and li arc the mass and moment of inertia, and _F!^ 
and Nf^ are the hydrodynamic force and torque exerted 
by the solvent on the beads, respectively [l^ [l^. 
and Gp are the random force and torque, respectively, 
due to thermal fiuctuations. The temperature of the sys- 
tem is defined such that the long-time diffusive motion of 
dispersed particles reproduces correct behavior [13, [H] . 

F^ represents the potential force due to direct inter- 
bead interactions, such as Coulombic and Lennard- Jones 
potentials. We use a bead-spring model as a model of 
a polymeric chain with a truncated Lennard- Jones po- 
tential and a finitely extensible nonlinear elastic (FENE) 
potential. The truncated Lennard- Jones interaction is 
expressed in terms of [/lj: 




+ e (r,j < 2scr), 



where rij = \Ri — Rj\. The parameter e characterizes 
the strength of the interactions, and a represents the 
diameter of the beads. Consecutive beads on a chain 
are connected by a FENE potential of the form 



FENE 



(r) 



-fcci?nlnn 



r 

Rq 



(6) 



where r = \Ri+i — Ri\, kc = 30e/cr^, and Rq = 1.5a. 
Ff is the constraint force acting on the ith bead due to 
the bond-angle constraints that cause the chain to form 
a straight line, and it is used only for the rigid rod case. 
This is given by 



d 



N 



dRih 



*a), (7) 



^„ = (a - 2)i?i - (a - l)i?2 



(a = 3,--- ,iV),(8) 



where 5'^ = is the constraint condition to be satis- 
fied, fla is a Lagrange multiplier associated with the 
intramolecular forces of the constraints chosen such that 
the constraint condition ^I'q, = is satisfied at a time t+h, 
where h is the time increment of a single simulation step. 
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B. SIMULATION 

Numerical simulations have been performed in three 
dimensions with periodic boundary conditions. The lat- 
tice spacing A is taken to be the unit of length. The 
unit of time is given by pf where 77 = 1 and pi = 1. 

The system size is x Ly x ^ 32 x 16 x 64. The 
other parameters include: fT = 4,^ = 2,e — 1,77 = 1, 
A/,; = 47raV3, iV = 5 and /i = 0.067. 

In the presented simulations, the Navier- Stokes equa- 
tion is discretized with a de-aliased Fourier spectral 
scheme in space and with a second-order Runge-Kutta 
scheme in time. To follow bead motions, the position, 
velocity and angular velocity of the beads arc integrated 
with the Adams-Bashforth scheme. 

At t = 0, the chain aligns along the x-axis, which is the 
shear direction. The run-time of our simulations is about 
3520/7. The range of ksT is 5.0 x 10"'' < ksT < 1.0, 
and that of 7 is 1.0 x 10"^ < 7 < 4.0 x lO^^ 

Prom the symmetry of the system, to analyze the tum- 
bling motion of a chain, we only have to consider the pro- 
jected tumbling motion on the x-y plane. We introduce 
the chain orientation angle (p, which is the angle between 
the a:-axis and the projected end-to-end vector on the x-y 
plane. 

Evidence of periodic tumbling motion should appear 
in the PSD per unit time [l, Q. Therefore, to investigate 
the spectral properties of the orientation, we calculate 
the PSD using a fast Fourier transform. The PSD is 
expressed as 

PSD{uj) = I / ip{t)cxp{iujt)dt\. (9) 



III. RESULTS 




V/Vp 



FIG. 1: (Color online) The behavior of PSD as a function of 
u/up for rigid rods at ksT = 0.0005 and 7 = 0.001 (blue □), 
ksT = 0.006 and 7 = 0.002 (red Q), and ksT = 0.01 and 
7 = 0.008 (green A). 




v/Vp 



In Figs. [T] and [U the PSDs of the chain orientation 
angle ip show a peak at a specific frequency i^p, as shown 
in previous studies Furthermore, the PSD data ob- 

tained at different conditions of shear rate and tempera- 
ture tend to lie on a single master curve if a normalized 
frequency v/vp is used. This is true for both rigid rods 
and flexible chains. This result implies that the tumbling 
motion of chains is fully characterized by Vp. 

In Fig. [3l we find that Vp follows the law Vp a j for 
7 < 0.02 at knT = 5.0 x lO"-*, for both the rigid rod 
and the flexible chain cases. Meanwhile, Fig. 2] shows 
Vp (X 7°-^** for 7 < 0.02 at fceT = 1.0 for both. 

Figure [3] also shows that, for 7 > 0.02, lyp is lower 
than the frequency expected from the law Vp (x j. P. 
Bagchi and S. Balachandar have reported that, at a finite 
Reynolds number approximately equal to I, the sphere 
rotation frequency in a linear shear flow decreases at a 
much slower rate than 7/2 [20|. The Reynolds number, 
i?e, is given by Re = pfja'^N/rj. In our paper. Re is 
equal to 1.6 at 7 = 0.02. We only consider I'p in the re- 
gion of 7 < 0.02, so the effects of finite Reynolds numbers 



FIG. 2: (Color online) The behavior of PSD as a function of 
u/up for flexible chains at ksT = 0.0005 and 7 = 0.001 (blue 
□), ksT = 0.006 and 7 = 0.002 (red Q), and ksT = 0.01 and 
7 = 0.008 (green A). 



do not influence the results. 

We arrange data sets using the Peclet number to con- 
sider the effect of competition between shear and fluctu- 
ation. The Peclet number is the dimensionless number 
that relates the rate of shear flow to the rate of thermal 
fluctuation. In our work, the Peclet number, Pe, and 
dimensionless frequency, F, are expressed as 



Pe 



F 



kBT 



GTrrja^Vp 
kBT 



(10) 
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FIG. 3: (Color online) The behavior of i/p as a function of 
7 at kuT = 0.0005. Rigid rod (green closed circle), flexible 
chain (red open circle) . The blue solid line corresponds to the 
law Up oc -y. Error bars are derived from the half band width 
of PSD. 



FIG. 5: (Color online) The behavior of _F as a function of Pe 
for a rigid rod. The orange solid and blue dotted lines were 
calculated from a least-squares fit of the data points. The 
orange solid line corresponds to Pe"'^^. The blue dotted line 
corresponds to Pe. 
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FIG. 4; (Color online) The behavior of Vp as a function of 7 
at /cbT = 1.00. Rigid rod (green closed circle), flexible chain 
(red open circle). The blue solid line corresponds to the law 
Vp (X 7^''^. Error bars are derived from the half band width 
of PSD. 



FIG. 6: (Color online) The behavior of -F as a function of Pe 
for a flexible chain. The orange solid and blue dotted lines 
were calculated from a least-squares fit of the data points. 
The orange solid line corresponds to Pe"'^*. The blue dotted 
line corresponds to Pe. 



We plotted the behavior of as a function of Pe for 
the rigid rod and flexible chain cases in Figs. [S] and [51 F 
was found to depend only on Pe because the data sets 
have the same value of F with the same value of Pe, even 
when the shear rates and temperatures are different. In 
the rigid rod case, F cx Pe°-^^ for Pe < 106 and F oc Pe 
for Pe > 106. In the flexible chain case, F (x Pe°-^^ for 
Pe < 156 and F oc Pe for Pe > 156. The behaviors of 
F for the rigid rod case and the flexible chain case are 
roughly equal, although the values of the Peclet numbers 
are different at the boundary where the exponent of Pe 
changes. 



We define the critical value at which the exponent of 
Pe drastically changes from almost 2/3 to 1 as the critical 
Peclet number, Pe^. F oc Pe^l^ for Pe < Pec; otherwise, 
F oc Pe. When fluctuations dominate the system, the 
exponent is nearly equal to 2/3. On the other hand, when 
shear flow dominates, the exponent is exactly equal to 1. 

In order to understand the behavior of F in the limit 
of Pe — > 00, we examine the behavior of in the limit 
of /cbT — > 0. As shown in Figs. [5] and El the behavior 
of F in the region of Pec < Pe is roughly equal to the 
behavior of v^, in the limit of k^T — >■ 0. Therefore, we 
believe that the proportional relation F oc Pe can be true 
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FIG. 7: (Color online) The behavior of Vp as a function of 7 
at ksT = 0. Rigid rod (red +), flexible chain (green x). The 
data sets for the rigid rods correspond to the law Up = O.I787. 
The data sets for the flexible chains correspond to the law 
Vp = 0.f907. 



across the entire region of Pcc < Pe. 



IV. DISCUSSION 
A. COMPARISON WITH OTHER RESULTS 

In our work, we calculate the tumbling motion of a 
single chain for 0.0005 < k^T < 1.00 and 0.001 < 7 < 
0.04. As reported in Jeffrey's paper [l^l, which treated 
non-Brownian particles, i^p follows the law oc 7 at 
/cbT = 5.0 X 10"^, as shown in Fig. H In previous 
papers that treated Brownian particles, i^p follows 
the law Vp (x 7^/'^ at ksT = 1.0, as shown in Fig. |31 

The dimensionless frequency, F, depends only on Pe. 
This Pe dependence on F can be described by a power 
law, F cx Pe". The exponent a drastically changes from 
2/3 to 1 at Pec] a equals 2/3 for Pe < Pec, while a 
equals 1 for Pec < Pe. In the case of a rigid rod with 
N ~ 5, Pec « 106, and Pec ~ 156 in the flexible chain 
case. When fluctuations dominate the system, F follows 
the law of F cx Pe^^^. On the other hand, when shear 
flow dominates, F follows the law of F cx Pe. 

Gerashchenko and Steinberg Q claim that there are 
two dynamical regimes of polymer motion at Wi 3> 1, de- 
pending on the polymer extension R. When R <C Rmax, 
where Rmax is the maximum polymer extension, the tum- 
bling frequency v is constant and independent of Wi. On 
the other hand, when R ~ Rmax, v is proportional to 
Wi^l^. In Fig. 4, we did not observe Vp to be indepen- 
dent of 7. We consider the chain length to be too short 
in our work. The chain cannot keep the coil state along 
the shear direction for long times. The chain is fully 



stretched at short notice and is always rotated, although 
Pe is small. 

Previously, Szymczak and Cieplak [2l| discussed the 
conformational dynamics of a single long protein in shear 
flow and found two characteristic tumbling frequencies, 
vi and ■ They showed that the higher frequency vi fol- 
lows the law of vi cx 7; however, the lower frequency 1^2 
follows the law of 1/2 cx 7^/'^ . When the protein is tightly 
packed, it essentially shows a spherical rotation in shear 
flow. As a result, vi is proportional to 7. The lower fre- 
quency, 1^2, corresponds to the stretching-collapse cycle; 
hence, V2 is proportional to 7^/'^. Although our results 
are similar to theirs, the phenomena in our system are es- 
sentially different from those in their works because rigid 
rods cannot fold. 

Davoudi and Schumacher analyzed the stretching 
of polymers in a turbulent flow. It is known that the 
polymers undergo a coil-stretch transition at Wi « 1/2 
in this system. For Wi < 1/2, polymers are in the coiled 
state, and their size distribution is stationary. In con- 
trast, for Wi > 1/2, the polymers are in the stretched 
state. Their stretching carries on until their lengths reach 
the finite extensibility limit or until turbulence stops the 
growth of the polymers. They found the maximum Lya- 
punov exponent to be A ^ j^^^. However, Chertkov ct al. 
reported A ^ 7^/'^ in their work, where A~^ is expressed 
as the mean stretching time scale. Davoudi and Schu- 
macher claimed also that their study could not be com- 
pared with the analytic results of Chertkov et al. They 
defined the shear time scale as Tg ~ 7^^ and the fluctu- 
ation time scale as Tf = D~^, where D is the strength 
of Gaussian fluctuation. Chertkov analyzed the polymer 
dynamics in the region of <^ Tf . Davoudi analyzed the 
polymer dynamics in the region of < Ts, where is 
the Kolmogorov time. Notably, Davoudi's work studies a 
different regime of polymer stretching than the analytic 
model of Chertkov's work. Chertkov et al. studied in 
the shear-dominated regime, whereas Davoudi and Schu- 
macher studied in the turbulence-dominated regime. 

Our work analyzed polymer dynamics in the region of 
Tg <C Tf . We do not consider the effect of turbulent flow. 
We only considered the region with particle Reynolds 
number Re < 1.6, so the effects of flnite Reynolds num- 
bers do not influence the results. In Davoudi's work, 
by contrast, the region of particle Reynolds number in 
the stretched state is estimated to be Re > 15, where 
Re ~ vL/v, V is kinetic viscosity, v is the root- mean- 
square of the turbulent velocity fluctuation, and L is the 
mean length of polymer. Therefore, our study cannot be 
compared with Davoudi's work. 



B. THEORETICAL ANALYSIS 

The geometry of the chain in our paper is depicted 
in Fig. [51 In Jeffrey's work [l^, the angle ip of non- 
Brownian rigid rods with flnite aspect ratios is governed 
by the equation 
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FIG. 8: The geometry of the present simulations. 
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FIG. 9: (Color online) Schematic diagram of two areas and 
the angle ipt- The angle ipt is the boundary angle of the area 
dominated by the first term and the area dominated by the 
second and third terms on the right-hand side of Eq. (|17p . 



In the case of D,. <^ 7, the dynamics of the angle (p be- 
come decoupled from the angle 9 between the end-to-end 
vector and the x ~ y plane because the angle 9 is ap- 
proximately zero. We can then write down the following 
equation: 



sin^ tp 



r 



(12) 



where the aspect ratio is r = L/a and L is the length of 
the chain. 

In our work, we consider the equation that governs 
the angle </3 of a Brownian rigid rod with a finite aspect 
ratio. To consider the diffusion of thermal fluctuation, we 
introduce white noise C(t) into Eq. p^ . We can then 
write down the following equation for a Brownian rigid 
rod: 



dt 



Dr 



cos^ 9 



C, (13) 



(14) 



where Dr is the rotational diffusion constant. On the 
basis of the shell model [2^, [23| , the rotational diffusion 
constant Dr for a rigid rod is calculated as 



dt"^ 



-7 



sm ip — 'j- 



1 



2^rC (17) 



For short times, the mean square displacement of ip{t) in 
time t is written as 



< \p{t) ~ ipiO)\'^ >^ 4Drt {for Drt<.l). (18) 

In order to analyze the tumbling motion of a single 
chain, we only have to consider the dynamics in the 
stochastic area, where the effect of thermal fluctuation 
is not negligible. In the case of Dr <C 7, the time re- 
quired to pass through the stochastic area is sufficiently 
larger than the time required to pass through the other 
area. The tumbling motion can only be understood by 
considering the stochastic area. The stochastic area is 
dominated by the second and third terms on the right- 
hand side of Eq. (|17p . The remaining area is dominated 
by the first term on the right-hand side of Eq. pT|) . Fig- 
ure [5] shows the two areas and the angle ipt , which is the 
boundary angle between the two areas. We call the sec- 
ond and third terms on the right-hand side of Eq. (fT7|) 
the shear and the fluctuation terms, respectively. 



Dr = 



3{ln r + d{r))kBT 
nrjL^ ' 



d{r) 



-0.662 



0.917 0.05 



(15) 



(16) 



In the shell model ixientioned above, the contour of 
the macromolecules of arbitrary shape is represented by 
a shell composed of many identical small beads. The shell 
model can be adequately modeled by decreasing the size 
of the beads. 



1. INFINITE ASPECT RATIO CASE 

First, wc consider Eq. (|17|) in the limit where r — > 
00. This limit is consistent with treating a chain as an 
end-to-end vector. We can rewrite Eq. ([T7|) as (yi = 
— 7sin^(/3 -|- 2yyDrC, in which the shear term does not 
exist. This equation is identical to the equation used in 
previous works [^-B ■ 

Previous works [J-Q have reported pt ~ (Dr/j)^^^- 
In the region of p>t < ^ the shear flow rapidly 

rotates the orientation of the chain from n to pt in a 
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time ti ^ In the region oi < (p < (pt, the chain 

orientation almost ahgns along the shear direction. Be- 
cause the effect of shear becomes sufficiently small in 
this region, the effect of thermal fluctuations only con- 
tributes to rotate the chain orientation from ipt to in 
time tr = iff/4:Dr oc as calculated with Eq. (fT8| . 

In the case of tr > ti, tr dominates the chain tumbling 
time. We can consider Vp cx t~^. Therefore, we conclude 
from Eqs. ^ and (HD) that F oc Pe^/^. 



2. FINITE ASPECT RATIO CASE 

Next, we consider Eq. ([T7|) with a flnitc aspect ratio, 
r. We can expect F to be proportional to Pe^/"^ when 
the fluctuation term dominates in the region < ip < tpt, 
as this case agrees with the limit of r — !• od. Additionally, 
we can expect F to be proportional to Pe when the shear 
term dominates in the region < (p < ipt because this 
case corresponds to the non-Brownian rigid rod. Pcc is 
deflned as the Peclet number at which the dominating 
term in the stochastic area changes from the fluctuation 
term to the shear term with increasing shear rate. 

From Eq. (fTT)) . we can write down the corresponding 
Fokkcr-Planck equation as 



dt 



d 

^•^'^ -0, 



(19) 



1 



1 



sin 95 + 7- 



1 



Wr^\P{t,pW 



where P(i, p>) is the PDF of the angle p> and is the 
probability flow. Each term in the braces of Eq. ([20)1 
corresponds to a respective term on the right-hand side 
ofEq. (HZD. 

Next, we focus on the stationary PDF P&t(}p) = = 
0, </3). Because dPst{'P')/dt = 0, namely, dJ^/dp — 
from Eq. (fT9)) . is 1^9- independent and constant: = 
-7[(r2 - 1) sin^ ip^P^tiPp) + Pst(<y5p)]/(r2 + 1), where p^ 
is the angle at which Pgt {p}) has a peak. We surmise that 
the angle pt satisfles the equation given by 



. (r'^-l)sinVP,t(y) 
7 r2 -f 1 



(21) 



The angle pt is calculated by solving Eq. (|22)) . Psx(}Pt) 
is expanded in powers l\p = pt — p-p. When we neglect 
the second order of /\p and higher, we can relate p)t to 

p-p as 




(23) 



It should be noted that pt,Pp <C 1. We can estimate pt 
with the angle pp. 

Next, we attempt to calculate the analytical form of 
ipp. In the case of Dr/j ^ 1, the formal solution for 
Psti^p) is given by 



PM = Ci rfV^exp ^-^J{i;,p) ] , (24) 



/(V',^)=^-(l 



1 



) sin?/'cos('0 — 2(y9), (25) 



where Ci is determined from the normalization condition, 
/p Pst(p)dp = 1. In the limit of r — > 00, it is known that 
Pp = {D^/^Y/^ 0. When r is flnite, pp is different from 
{D-c/^Y^^ because the shear term influences Pst{p)- 

To estimate the effect of the shear term, we introduce 
the times ^31 where ^3 represent the times 

required to pass through the region < 1^ < {D^/'yY^^, 
depending only on the shear term and the fluctuation 
term, respectively. <2 is estimated by dividing {D^/jY^^ 
by the shear term as 



(r^ + VfDi 

. 4 

73 



(26) 



and t^ is estimated by Eq. p8| as 



*3 = 



(27) 



Because the Brownian chain is rotated by both the 
shear term and the fluctuation terms, we expect that p>p 
is given by 



Pp^C(r,Pe)[^ 

1 



(28) 



By substituting the expressed by pp into Eq. (PT|). 
the relation of pt to p^p is given as 



. (r^ - 1) sin^ ppP^t{pp) + PstjPp) 

^. (r2-l)sin' PtPA^t) 
r-^ + 1 



C (r, Pe) 



1 + 



18(ln r+£i(r)) 



Pe 



(29) 



For large 7/-Dr where the present theoretical analysis 
is valid, C {r,Pe) ^ ^2/^3 finally approaches zero. We 
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FIG. 10: (Color online) The behavior of ipt as a function of 
Pe. ifi^i is ipt calculated with ipp, obtained from numerical 
integration (red +). The green solid line corresponds to the 
behavior of Lp^n , where Lp^n is (pt, obtained from ipp = [42/(^2 + 

confirmed, however, that ^2/^3 still remains finite (^ 1) 
around Pcc - 100, where j/Dr « 1000. 

Figure [TO] shows the behaviors of (p^' and ip^" as a func- 
tion of Pe, where ip^i is ift calculated with Eq. ([23]) and 
ipp, which is obtained from the numerical integration of 
Eq. ([24]), and (pf." is 1^94 calculated with Eq. ([23]) and (fp, 
which is expressed as Eq. (j28p . In this figure, it is shown 
that Eq. ([28]) is established because the behavior of ip^' 
agrees well with ip^" for Pe < 100. We can obtain ipt 
by substituting Eq. ([21]) into Eq. Therefore, (pt is 

given by 



ip, = C{r,Pe)(-^Y , (30) 



We consider that F is proportional to Pe^/^ when the 
effect of thermal fluctuation is more significant than the 
effect of shear flow in the region of < </? < and that 
F is proportional to Pe in the opposite case. t2 is the 
time required to pass through the region Q < ip < Lpt hy 
the shear term, 

= (li±i)£W (32, 

ts is the time required to pass through the region < 



<y9 < (y3t by thc fluctuation term, 

h = (33) 
41) ■'7 3 

It is thought that Pcc is thc Peclet number that satis- 
fies t2 = ts- If wc obtain thc value of r, we know the value 
of Pec because both t2 and are functions of Pe and r. 
Our results show that Pec « 115 at r = 5. In thc nu- 
merical results obtained from our work. Pec = 106. The 
analytical result agrees well with our numerical result. 
Moreover, the numerical condition Dr ^ 7 is satisfied 
because D^/j ss 10"^ < 1 at Pe sa 100. Therefore, the 
considerations in this section are reasonable in thc region 
near Pec- 

The considerations in this section a gree well with those 
of previous experimental results 0-11, 113, [111 ■ In experi- 
mental works that measured the frequencies of DNA rota- 
tion [i-[l|l3|, the DNA molecules contained roughly 400 
persistence lengths. The persistence length is thought to 
correspond to r. Thus, we conclude that Pec « 516. 
These experiments were carried out in the region of 
Pe <C Pec, and F is proportional to Pe^^^. In exper- 
imental work that measured the frequencies of freely ro- 
tating rigid dumbbells [3] ,the aspect ratio r of thc rigid 
dumbbell corresponds to 2 and Pec ~ 60. The experi- 
ments were carried out in the region of Pec <C Pe, and 
F was proportional to Pe. From these results, the con- 
siderations in this section are reasonable. 



V. CONCLUSION 

In our work, we calculated the tumbling motion of 
a single chain using an SPM that takes into account 
thermal fluctuations and hydrodynamic interactions for 
0.0005 < k^T < 1.00 and 0.001 < 7 < 0.02. We conclude 
that the dimensionless frequency, F , depends only on Pe. 
Thc dependence of F can be described by a power law 
F (X Pe". The exponent a sharply changes from 2/3 to 1 
on Pec. In the case of a rigid rod with TV = 5, Pec ~ 106, 
and in the case of a flexible chain with iV = 5, Pec ~ 156. 
The behavior of F for both cases is similar, while only 
the values of Pec are different from each other. 

We have presented F to be proportional to Pe^^^ when 
thc third term on the right-hand side of Eq. ([T7]) dom- 
inates in the region < < ipt, and F is proportional 
to Pe when the second term of Eq. ([T7]) dominates in 
the region < (p < (pt. We have estimated the angle ipt 
at which the fist term of J^, expressed as Eq. ([20]) . is 
comparable to thc sum of the second and third terms of 
Jip, expressed as Eq. ([^D]) . 

A proposed mechanism for this exponent change is that 
the effect of thermal fluctuation is more significant than 
the effect of shear flow only for < ip < ipt, whereas in 
the other case, the effect of thermal fluctuation is negli- 
gible. The former contribution leads to F cx Pe^/^, and 
thc latter contribution leads to F cx Pe. 
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